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The critical exponent /3 = 0.16 ± 0.02 for the random-field Ising model order parameter is de- 
termined using extinction-free magnetic x-ray scattering for Feo.85Zno.15F2 in magnetic fields of 10 
and 11 T. The observed value is consistent with other experimental random-field critical exponents, 
but disagrees sharply with Monte Carlo and exact ground state calculations on finite-sized systems. 

PACS numbers: 61.10.Nz, 75.40.Cx, 75.50.Ee, 75.50.Lk 



The Ising model is perhaps the most important model 
in statistical physics, with many applications that go well 
beyond the realm of physics. In 1944, the pure two- 
dimensional (d = 2) model was solved exactly by On- 
sager 0. For the pure d — 3 Ising model, although 
various calculations and computer simulation techniques 
have proven extremely useful, no exact results are avail- 
able. Nevertheless, there is extremely good agreement 
among theory, simulations, and experiments Hence, 
the model can be considered well understood. Through 
the years, the intellectual value of the Ising model has 
grown, particularly as a model of disorder. One of the 
most important of these models of disorder occurs when 
a random field is imposed which couples directly to the 
order parameter of the system. The most studied re- 
alization of this random-field Ising model (RFIM) is the 
diluted antiferromagnet in an applied magnetic field. Un- 
like the pure d — 3 Ising model, there is, so far, poor 
agreement between theory and simulations on the one 
hand and experiments on the other Q|. This has moti- 
vated us to measure the critical behavior of the staggered 
magnetization, the antiferromagnetic order parameter 



M s = Mot 



(1) 



where t = (T C {H) - T) /T C (H), of the dilute Ising antifer- 
romagnet in an applied field, since this is one of the most 
valuable yet also least characterized aspects of the exper- 
imental system. Our result for the order parameter expo- 
nent j3 provides an important quantitative experimental 
contribution toward a comprehensive understanding of 
the RFIM. 

The dilute, insulating antiferromagnet Fe x Zni_ x F2 has 
a large single-ion anisotropy and is an extensively studied 
§ d=3 RFIM realization f|, ||- Nevertheless, prior at- 
tempts to determine the critical behavior of the order pa- 
rameter have been unsuccessful. This may be surprising 
since, in principle, one only needs to measure the temper- 
ature dependence of the Bragg scattering intensity, I B , 



which is proportional to (M s ) 2 , with a magnetic field H 
applied along the c-axis, the spin ordering direction. For 
two reasons, such measurements have proven very diffi- 
cult in practice. First, neutron scattering on these high- 
quality crystals suffers greatly from the effects of extinc- 
tion; the beam, upon transmission through the crystal, 
is depleted of neutrons satisfying the Bragg condition, 
resulting in the saturation of the measured value of Ib ■ 
Second, for magnetic concentrations x below the vacancy 
percolation threshold [|, |(|, x v — 0.754, domain forma- 
tion obscures the RFIM critical behavior below the tran- 
sition at T C (H) in Fe x Zn!_ x F 2 and its less anisotropic 
isomorph Mn x Zni_ x F2. Although the domains may be 
internally well ordered, Ib will be greatly diminished if 
the characteristic length scale for the domain structure is 
smaller than that of the spectrometer resolution; we will 
refer to this as micro-domain structure and it has been 
studied extensively in previous works Q|. Under severe 
extinction conditions, domain structure may relieve ex- 
tinction and actually cause Ib to increase. Whether do- 
main structure forms or not, the Bragg scattering cross 
section will be decreased by thermal disorder as the tran- 
sition is approached. In x-ray scattering, since the mag- 
netic scattering cross section is relatively small, the scat- 
tering intensity, obtained in a reflection geometry, does 
not suffer from extinction, as extensively discussed pre- 
viously |^|, ^] . The use of extinction- free magnetic x-ray 
scattering, and of a crystal with x > x v to avoid micro- 
domains, has allowed us to accurately characterize the 
order parameter critical behavior in Feo.85Zno.15F2. 

The magnetic x-ray scattering technique was employed 
for MnF 2 for H — by Goldman et al. § , and was then 
applied to Mn x Zm_ x F 2 with H > by Hill et al. [§. 
Whereas the H = study yielded the exponent (5 con- 
sistent with the d = 3 Ising model, the latter did not 
reveal the universal RFIM behavior ||, which would be 
consistent with x < x v ; the H = transition tempera- 
ture Jl(J and its field dependence are consistent with 
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the concentration being very close to or below x v . That 
the transition in Ref. |)j was obscured by micro-domains 
H is suggested by the zero slope of M s 2 versus T as 
T — > T C (H). The present magnetic x-ray scattering mea- 
surements use Feo.85Zno.15F2, for which x is well above 

Xy . 

The measurements were made at the new high-field 
magnet facility on beam line 7-2 of the Stanford Syn- 
chrotron Radiation Laboratory. A monochromatic x- 
ray beam was obtained from the wiggler spectrum via 
a Si(lll) double-crystal monochromator. X-ray energies 
between 14 and 13.5 keV were used, which resulted in 
a penetration depth of about 60 fim. The energy was 
tuned to minimize energy-sensitive multiple scattering 
P). The sample had a finely polished face, a few mm 2 
in area, with the a-axis perpendicular to the polished 
face and the c-axis along the vertical field. The tem- 
perature of the crystal, mounted in a He atmosphere, 
was stable to approximately 10 mK. The transition tem- 
perature for H = was measured to be Tn — 66.7 K, 
consistent with birefringence measurements on the same 
sample 111] and with a concentration x = 0.85 |l(J. It 
has been shown that the antiferromagnetic transition at 
this magnetic concentration is stable at fields as high as 
H = 18 T g. For H = 10 and 11 T, the transitions 
are at T c — 64.2 and 64.0 K, respectively. The lattice 
constants of the sample were determined to be approx- 
imately a = 4.68 A and c = 3.27 A near the transition 
temperature. The half-widths-at-half-maximum for the 
Bragg peaks were 4 x 10~ 4 , 4 x 10 -3 , and 4 x 10 -3 recipro- 
cal lattice units for the transverse, longitudinal and verti- 
cal directions, respectively, at the (100) magnetic Bragg 
point, about which transverse H = and H = 11 T 
and longitudinal H = 10 T scans were obtained. The 
sample was remounted between measurements at differ- 
ent fields, and therefore we normalized intensities using 
scans at T = 47 K. Three conventional thermal-cycling 
procedures were employed. In ZFC, the sample is cooled 
in zero field below T C (H) : the field is applied, and the 
sample is warmed through T C (H), waiting at each tem- 
perature at least 20 min before taking data to let the 
temperature and system stabilize. In FC, the sample 
is cooled through T C (H) in the field, taking data as in 
ZFC. Field- heated (FH) data were taken by heating in 
the field after FC. The scans typically consisted of 41 
points, about 15 of which covered the Bragg peak. At 
each point, the intensity was counted for 30 to 45 sec- 
onds, depending on the temperature of the scan. 

Figure 1 shows the Bragg intensity for H = and 10 T 
versus temperature, with the q and T independent back- 
grounds subtracted, where q is the distance in reciprocal 
space from the (100) antiferromagnetic Bragg point. The 
background depends on the precise experimental config- 
uration, but not on the thermal cycling used to collect 
data, and is mostly from sources other than the crystal 
itself. For comparison of the background to the Bragg 
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FIG. 1: The ZFC Bragg intensity, Ib, in arbitrary units (a.u.), 
versus T for H = and H = 10 T, with the T-independent 
background intensity subtracted. The square symbols in in- 
sets a) and b) show the critical scattering contributions to the 
x-ray intensities for H — and H = 10 T, respectively, de- 
termined from neutron scattering measurements as described 
in the text. 

signal, typical background counts for the H = 11 T scans 
were eight counts per second whereas the q — intensity 
was 160 counts per second at T = 47 K. Above the transi- 
tion, the scattering intensity results only from the critical 
scattering and goes to zero well above T C (H), indicating 
that there are no contributions from multiple scattering 
to the Bragg intensity. To determine the critical scat- 
tering for H = 0, 10 and 11 T, neutron scattering line 
shapes, obtained with a sample of nearly the same mag- 
netic concentration 1T3] using a previously described pro- 
cedure H] , were folded with the x-ray resolution, and the 
overall q = amplitude was adjusted to fit the H = 10 T 
data above T C (H). Insets a) and b) in Fig. 1 show the 
critical scattering contributions for H = and 10 T, re- 
spectively. As a result of the high momentum resolution 
of the x-ray technique, the critical scattering contribu- 
tions, which are nearly Lorentzian for H = |l5|] , are 
almost negligible (Fig. la)). For H > 0, however, the 
critical scattering has a much larger q dependence Jl3] 
at small q. Consequently, a small contribution to the 
q = scattering is more discernible for the H = 10 T 
data. These contributions were subtracted from all the 
data before determining the order parameter exponent. 

Although neutron scattering measurements using 
Feo.g5Zno.15F2 |H) and Feo.93Zno.07F2 @ show no evi- 
dence for micro-domain formation in the critical scatter- 
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ing, H > hysteresis in Ib is evident, with the FC in- 
tensities larger than the ZFC ones, a result of extinction. 
The x-ray Bragg scattering also shows hysteresis, but in 
this extinction-free case the ZFC data are higher in inten- 
sity. In the random field region, the FC intensity is fac- 
tor of 4.1 smaller than ZFC intensity. This ratio depends 
slightly on the cooling rate used in obtaining the FC data. 
The ZFC data are rate independent. For H = 11 T, the 
corresponding ratio is approximately 4.0. FH data were 
intermediate between the ZFC and FC curves. We note 
that specific heat measurements also show hysteresis very 
close to T C (H) at this concentration |lj, |l6). Such hys- 
teresis is likely a result of the extremely slow activated 
RFIM dynamics and possibly represents the fact that FC 
long-range ordering must take place while traversing the 
transition, where activated dynamics plays the greatest 
role. The logarithmically slow relaxation associated with 
activated dynamics |l7], [l8] severely limits the ability of 
the system to equilibrate extremely close to T C (H). This 
limits the formation of long-range order upon FC. ZFC 
data, on the other hand, are obtained without approach- 
ing T C (H) except when the order parameter is already 
very small, and thus do not visibly suffer from the slow 
dynamics. Moreover, no dependence on the rate of tem- 
perature change was observed upon ZFC. Hence, we be- 
lieve the ZFC data represent the correct order parameter 
measurement. Various measurements near the transition 
at this concentration have yielded critical behavior in- 
dicative of a second-order phase transition. There is also 
no measureable latent-heat in specific heat critical behav- 
ior measurements fL2|| . Therefore, this appears to be a 
second-order transition, although an extraordinary one. 

The normalized Bragg intensity curves in Fig. 1 clearly 
approach T C {H) vertically. This is characteristic of ex- 
periments fl|, |lj, [l5| and simulations Hi for x > x v 

and 
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and in stark contrast with experiments 
simulations Q for x < x v , where Ib approaches T C (H) 
horizontally. The latter behavior is attributable to micro- 
domain formation, which is energetically favorable when 
the vacancies percolate through the crystal, as shown in 
Monte Carlo simulations ||. 

Figure 2 shows the logarithm of Ib, with the con- 
stant background and critical scattering contributions 
subtracted, for H = 0, 10, and 11 T, versus the loga- 
rithm of t. The values of T C {H) were determined from 
fits to the data. For 0.0007 < t < 0.03 and H = 0, we 
find (3 = 0.35 ± 0.02 (solid line), which agrees well with 
several experimental and theoretical determinations for 
the random-exchange Ising model B. For H = 10 and 
11 T, a crossover from random-exchange to RFIM critical 
behavior occurs near t — 0.03, consistent with birefrin- 
gence measurements ||l2| , and the data can be fit to a sin- 
gle power law only in the range 0.0001 < t < 0.03. The 
fits over this range yield the exponent j3 = 0.16 ±0.02 for 
the combined H = 10 and 11 T data and are indicated by 
the parallel solid lines in Fig. 2. A less sophisticated data 
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FIG. 2: The same ZFC data as in Fig. 1 as well as data taken 
at H — 11 T, corrected for the critical scattering contribution, 
plotted as the logarithm of the intensity versus the logarithm 
of t. The solid lines for H = 10 and 11 T indicate RFIM 
behavior with /3 = 0.16, while the solid line for H — reflects 
conventional random-exchange behavior (/3 = 0.35). 



analysis, which did not correct for the critical scattering 
contribution, resulted in a value for the order parame- 
ter exponent that is larger by 0.02, still within the error 
bars. The correction should, of course, be done in order 
to obtain the correct value of (3. The slope of the data 
at large t changes with H since the definition of the re- 
duced temperature involves T C (H) for H > and not the 
H = transition temperature Tjv. 

Through the Rushbrooke scaling relation 



2/3 + 7 + a > 2 



(2) 



which is usually satisfied as an equality, (3 is related to the 
universal critical exponents a (for the specific heat) and 7 
(for the staggered susceptibility) of the d = 3 RFIM. The 
experimentally determined specific heat peak is nearly 
logarithmic and very symmetric close to T C {H), consis- 
tent with a w (l^, fl6]|. Neutron scattering analyses 
p3| , fu] yield values in the 1.45 < 7 < 1.65 range. There- 
fore, the experimental value f3 f=a 0.16 is fairly consistent 
with Rushbrooke scaling, taking the upper limit of 7 and 
a = 0. 

A very recent NMR study J2(| of the order parame- 
ter in the effective short-range interaction random-field 
ferroelectric Sr .6i-xCe x Ba .39Nb2O 6 , with x — 0.0066, 
yielded (3 = 0.14±0.03, consistent with the present result 
for Feo.85Zng.15F2. It was obtained, however, in the pres- 
ence of micro-domain structure. Apparently, the more 
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local NMR probe is less sensitive than the Bragg scatter- 
ing techniques to the formation of micro-domains. This 
suggests that domain formation does not preclude a fairly 
sharp RFIM-like phase transition and a measurement of 
its order parameter , but only prevents measurements 
of the order parameter through scattering experiments 
for x < x v . 

There does not exist a set of theoretical results that 
are consistent with all the experiments Monte Carlo 
p|| and exact ground state calculations [^3| yield very 
small values for [3 and large, negative values for a. 
Other numerical and scaling analyses p4| yield a close to 
zero, consistent with experiments fl2|, |16||, but also yield 
v = 1.37 ± 0.09, much larger than experimental value 
v = 1.05 ± 0.01 (R§ 0. Another recent work || yields 
a and (3 close to zero. One Monte Carlo study, on a large 
lattice and with less assurance of equilibrium than other 
simulations, yielded f3 = 0.25 ±0.03 ^6|. Since con- 
sistency among numerical and experimental exponents 
continues to elude us, a comprehensive understanding of 
the d = 3 RFIM is yet to be achieved. The determina- 
tion of the order parameter exponent presented here is 
an important quantitative contribution in this direction. 
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